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Abstract 


The gravitational influence of a planet on a nearby disk provides a powerful tool for detecting and 
studying extrasolar planetary systems. Here we demonstrate that gaps can be opened in dynamically 
cold debris disks at the mean-motion resonances of an orbiting planet. The gaps are opened away 
from the orbit of the planet itself, revealing that not all disk gaps need contain a planetary body. 
These gaps are large and deep enough to be detectable in resolved disk images for a wide range of 
reasonable disk-planet parameters, though we are not aware of any such gaps detected to date. The 
gap shape and size are diagnostic of the planet location, eccentricity and mass, and allow one to 
infer the existence of unseen planets, as well as many important parameters of both seen and unseen 
planets in these systems. We present expressions to allow the planetary mass and semimajor axis to 
be calculated from observed gap width and location. 

Subject headings: planet-disk interactions - planets and satellites: detection 


1. INTRODUCTION 

The detection of excess infrared radiation from a num¬ 
ber of stars indicate the presence of debris disks around 
them, some of which, such as those around Fomalhaut 
and Beta Pictoris, have been observationally resolved. 
The first extrasolar debris disk ever discovered was found 
around Vega using IRAS\ from the thermal emission of 
circumstellar dust which revea led a strong infrared excess 
beyond 12 fim (Aumairn et al, 19841) . The excess radia¬ 
tion was immediately linked to the possible presence of 
solid dust particles with radii greater than 0.12 cm and 
temperatures of approximately 85 K, located at a mean 
distance of 85 AU from the nearby main sequence star 
a Lyrae and which were believed to be debris left-overs 
from the formation of this stellar system jAumann et al.l 
[1984) . This discovery sparked interest in studying de¬ 
bris disks outside our own solar system and the first 
optical image of an exosolar debris disk emerged later 
in the same year. Using ground-based optical coronag- 
raphy, Smith & Terrile were able to directly image a 
flattened disk of cold, solid material around ft Pictoris 

([Smith h Terrildfl984|) . 

Interactions between planets and disks have been stud¬ 
ied in great detail for gas-rich disks with the main mo¬ 
tivation being u nderstanding planet formation (see the 
review paper by lKley fc Nelson|l2Q12|) . The study of ge¬ 
ometric structures in protoplanetary disks has offered 
insights into the formation and evolution of planetary 
systems. Such structures, which are mostly believed to 
be signposts of planet formation, reveal themselves as 
density variations across the disk and have been inferred 
from Spectral Energy Distributions and high-resolution 
millimeter and submillimeter interferometry observa¬ 
tions. There is strong theoretical and observational ev¬ 


idence for gaps (see for instance, iDebes et~aLl I2013f) as 
well as density enhancements that appear as complex 
features such as spiral patterns (e.g . Brown et alj 12009 ; 
Vluto et al.||2012f [Juhasz ct al.|20l4) and dust traps (e.g. 
Isella et all 12013: (van Her Marel et aLl 120131) in proto- 
planetary disks and that can create asymmetric struc¬ 
tures. Observations of protoplanetary disks have re¬ 
vealed that asymmetric disks are common. Such asym¬ 
metries are interpreted as being either due to density 
perturbations of a stellar or planetary companion (e.g. 
IKraus ctld] |2013f) or having a ge ometric nature (e.g. 


Brown etaD 120091; [Jang-Condcll & Turned [2013) . Inves¬ 
tigating the sources of density enhancements and deple¬ 
tions can result in a better understanding of the processes 
of formation and evolution of single and multiple plane¬ 
tary systems in gas-rich disks. 

In a protoplanetary disk, drag against the gas causes 
solid particles to collapse into a dynamically cold disk in 
which the particle orbital eccentricities and inclinations 
are very low. The gas is eventually blown away by stel¬ 
lar radiation once the star is born and leaves behind a 


nearly circular, coplanar solid particle disk.This disk may 
include planets as well as smaller planetesimals or other 
bodies, like the solar system’s asteroid belt. Dynami¬ 
cal interactions of planets with these second generation 
disks have not been studied as extensively as gas- rich 
disks. In this paper, we shall focus on this later stage in 
which little or no gas remains. Su ch gas-poor disks in- 
clude for example, the F omalhaut (ICataldi et alj [20151), 

S dWilner et alJI2002f) and (3 Pictoris (jKalas & Jewitti 
disks. Structures in the Fomalhaut disk, for exam¬ 
ple, cannot be due t o gas-dominated pro cesses because of 
its low gas content ([Cataldi et alJl2Q15Q . but must arise 
from other processes. 

We wall show that the interaction of a planet with a 
debris disk can create structures that are not radially 
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symmetric about the star. Non-axisymmetric structures 
have been commonly observed in debris disks. For in¬ 
stance, observations of the debris disk around 3 Pictoris 
reve aled a warp in the inner disk around ~ 7 0 — 150 AU 
(e.g. Burrow s et al. Il995i iHeap et alj l20QQf) . Dynami¬ 
cal modeling had suggested that t he warp could be ex¬ 
plained by a misalign ed planet (e.g. IMouillet et al]|19971; 
lAugercau et al.l 12001ft which was later confirmed when 
a 9 ± 3 Jupiter-mass planet was found on an incline d 
orbit 8-9 AU from the central star (|Lagrange et al.l 
120101) . Asymmetries in debris disks have also been at¬ 
tributed to mean-motion resonances (MMRs). For in¬ 
sta nce, A-body simula t ions o f a collisional debris disk 
by iNesvold h Kuchncr^ (]2015l) show a peak in the disk’s 
surface brightness at 1:1 MMR. with the planet. A sec¬ 
ond peak is observed in their simulation of a 3 Jupiter- 
mass planet and falls between the 3:2 and 2:1 MMR 
with the planet, indicating a depletion of planetesi- 
inals at the two resonances. Formation of overdensi¬ 
ties in debris d isks has also been explained by migrat¬ 
ing dust (e.g. Wilner et al.l 120025 or migrating plan¬ 
ets (e.g. [Wyattj 12003ft . While migrating inward due to 
Poynting-Robertson (PR) drag, dust grains can become 
trapped in resonances with a planet interior to their or¬ 
bit or get scattered, thus forming asymmetric structures 
in the disk. The same scenario applies to particles that 
are captured in a migrating planet's resonances. 

Structures in dust disks in which particles are strongly 
affected by the radiation pressure of the central star as 
well as the PR. drag h ave been investigated by some au¬ 
thors (see for instance iK uchner &: Holmanl 12003 ; IWvattl 
|2006t IKrivov et alj 1200 7). However, not much empha¬ 
sis has been placed on the dynamical interactions of 
planets with planetesimal belts with regards to under¬ 
standing of the MMR gaps and how they can be used 
to extract info rmation abou t the p lanets causing them. 
For instance, IChiang et al.l d2009|) briefly discuss gap 
formation at MMRs with a planet interior to the Fo- 
malhaut disk in an attempt to constrain the mass of 
Fom b, but do not provide further details on how such 
gaps could yield measurements of planetary parameters. 
On the other han d, A-body numerical simulations by 
IReche ct"al] (1200815 show gap structures similar to what 
we shall discuss in the present work, but were not ad¬ 
dressed by the authors. Here we address the question of 
what structure might be induced in a gas-poor extraso¬ 
lar planetesimal disk by a non-migrating planet, and in 
particular what observational signatures might indicate 
the presence of an unseen planet in the system and how 
they can be used to constrain the planet’s mass and or¬ 
bital parameters. It must be noted that for the rest of 
this paper, whenever the term “debris disk” is used, it 
means disks with particles having size distributions in the 
range of about 1 m to 100 km that are nearly unaffected 
by the central star’s radiation pressure. Such particles 
can be gravitationally perturbed due to the presence of 
one or more planets in the system which could result in 
either their removal (like the Kirkwood Gaps in our solar 
system) or accumulations (like the Hilda family of main 
belt asteroids). Studying these resonances is not only 
an indirect way of detecting unseen extrasolar planets, 
it can also help put constraints on some parameters of 
the planets which are creating them, such as their mass, 
semimajor axis and eccentricity. 


The astonishingly detailed image of a disk around HL 
Tau, a Sun-like star approximately 450 light-years away 
in the constellation Taurus that was recently obtained by 
the Atacama Large Millimeter Array (ALMA), is a per¬ 
fect exa mple of structures being formed in disks due to 
planets (NRAO l2Q14ft . Although the disk imaged around 
HL Tau is a protoplanetary disk different from the second 
generation debris disks we concentrate on in the present 
work, the ALMA image of HL Tau illustrates the increas¬ 
ing resolving power that can be achieved with state of the 
art telescopes. With ALMA soon starting its full oper¬ 
ation, more detailed images of debris disks will become 
available, revealing more and more detailed structures 
such as those discussed here. 

This paper is organized as follows: In Section [2J we 
briefly describe the dynamics of MMRs and the theoreti¬ 
cal calculations of their maximum libration widths. Our 
method is presented in Section®in which we explain the 
code that is used to generate the simulations as well as 
the initial conditions. We present the different results 
obtained for ulterior versus exterior MMRs in Sectional 
followed by discussions in Section® Finally, a summary 
and conclusions are provided in Section® 

2. THE DYNAMICS OF MEAN-MOTION RESONANCES 

The existence of gaps in the solar system’s asteroid belt 
was first noted by American astronomer Daniel Kirk¬ 
wood in 1867 who saw non-uniformities in the number 
distribution of asteroids in the main belt as a function of 
semima ior axis, w ith some ranges having few or no as¬ 
teroids ([Kirkwoodl 1867ft . We now understand that grav¬ 
itational interactions between the asteroids and Jupiter 
result in the removed of planetesimals from these orbits, 
making gaps to appear in the disk where a particle's or¬ 
bital period would be a simple fraction of that of Jupiter’s 
([Murray fe Dermott] 119995 . The close link between or¬ 
bital period of an asteroid and its semimajor axis means 
that these “MMRs” occur over narrow ranges of semi¬ 
major axes, often depleting them of their original com¬ 
plement of bodies. 

However, if one were to take an image of our asteroid 
belt from outside the solar system, the Kirkwood gaps 
would not be observable due to the eccentricity of the 
asteroids blurring the edges of the gaps (see Figure 0. 
Though asteroids are removed from a number of reso¬ 
nances in the main belt, the eccentricities of the remain¬ 
ing asteroids are large enough to blur the edges of the 
resonances and make the gaps invisible. We will show 
that this is not always the case, and that under realistic 
conditions, visible gaps can be opened in particle disks by 
resonances. This means that gaps in extrasolar disks do 
not necessarily contain planets, as often assumed. Plan¬ 
ets can also generate other resonant structures that both 
indicate their presence and provide diagnostic informa¬ 
tion about their mass, eccentricity and position. There¬ 
fore, not only is the observation of gaps in debris disks an 
indirect way of detecting undiscovered extrasolar planets, 
it can also be used to constrain some parameters of the 
perturbing body if the gap widths and locations can be 
measured. 

Two objects are said to be in MMR. if the following 
relation holds: 

n p + q 
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Figure 1. O rbits of known solar system asteroids interior to Jupiter's as of 2015 June 27 plotted using data obtained from the Asteroids 
Dynamic Site (AstDyS120ir>f) . No Kirkwood gap due to MMR with Jupiter (shown with the black circle) can be seen in the disk (top panel) 
or the histogram (bottom panel) which shows the distribution of asteroids per distance from the Sim. This is due to the eccentricities of the 
asteroids that bring them in and out of the gaps. The theoretical locations of the three strongest resonances are shown on the histogram 
with vertical dotted lines. 


where n and n! are their mean-motions (= with T 
the orbital period), and p and q are positive integers, 
with q denoting the order of resonance. The primed and 
unprimed quantities are the orbital elements of the parti¬ 
cle being perturbed (the “asteroid”) and the perturbing 
body (“the planet”), respectively. Equation[T]can, equiv¬ 
alently, be written in terms of the two objects’ semimajor 
axes, a and a', as: 



where k = +1 for exterior resonance (i.e. a' > a) and 
k. = — 1 for interior resonance (i.e. a' < a). 

Following the dis c ussion and derivations presented in 
1 Murray fc Dermott] dl999() . to lowest order the resonant 
argument of the disturbing function, <f>, can be written 
as: 


where A and X' are the mean longitudes and cj and w’ 
are the arguments of periapse. Also, j\ = p + q , j% = —p 
and j3 and j± are either zero or — q , depending on the 
relative locations of the two objects. 

When a particle in a debris disk orbits in a MMR with 
a perturbing body, such as a planet, its orbit is per¬ 
turbed in a consistent manner when the relative planet- 
asteroid geometry repeats itself. This often destabilizes 
the smaller body so that it either collides with the planet, 
crashes into the star, goes into a highly elliptical orbit or 
gets ejected. In either case, a gap forms in the disk whose 
width, Sa ' max , can be approximated by the maximum li- 
bration width of its resonance. This can be calculated 
analytically at low eccentricities using Equation 0) 
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where: 

«-(s^) *«">, W 

with G denoting the Universal Gravitational Constant, 
M the mass of the perturbing body, fd(&) a term 
containing Laplace coefficients and coming from three- 
body expansion of the two objects’ orbital elements and 
a = Equa tion U is given as Eq. (8.76) in 

I Murray fe DermottJ ([19990 . Figure shows the maxi¬ 
mum libration widths at various eccentricities for the 2:1, 
3:2, and 3:1 interior resonances with Jupiter. It is clear 
from this figure that in the case of first order resonances, 
the libration widths increase with eccentricity, except at 
very low eccentricities for which they become very large. 

A MMR affects a narrow range of semimajor axes. The 
resulting resonant gaps will be visible in an optical or 
other telescopic image of the disk if (1) the resonance 
tends to destabilize particles (often but not always the 
case), (2) the eccentricity of the particles near the edges 
of the resonance is small enough to keep them from blur¬ 
ring the edges of the gap, and (3) the width of the gap 
is not smaller than the resolution of the telescope. The 
radial excursion, Ar, of the particles in the disk, given by 
Ar ~ 2ae, should be smaller than the resonance width. 
Thus the formation of MMR gaps is favored in dynami¬ 
cally cold particle disks, which have small eccentricities. 
It is the relatively high average eccentricities of asteroids 
that obscure the Kirkwood gaps in our own main asteroid 
belt. From Figure [2] it is also clear that very low eccen¬ 
tricity disks may have very large resonant gaps since the 
resonance width increases sharply as the disk average ec¬ 
centricity, e' avc , approaches zero; hence very dynamically 
cold disks are good candidates for observing particularly 
large resonant gaps. 

3. SIMULATIONS 

3.1. The Method 

Our simulations are performed with a symplectic 
integrator based on the Wisdora-Holman algorithm 
dWisdom fc Holman||1991[) . A fixed timestep of 50 days 
is used for all simulations. The output is recorded at 
10,000-year intervals and the total simulation time is 
taken to be 1 million years in length unless otherwise 
noted. A single planet on a circular or slightly elliptical 
orbit perturbs the disk, which orbits a 1 solar-mass star. 
Particles are removed if they have a close encounter with 
the planet or reach a distance less than 10 solar radii or 
greater than 1000 AU. 

The gravitational effects of the star and planet are in¬ 
cluded but interactions between the particles themselves 
are ignored. Thus our simulations are applicable to low- 
mass debris disks, where the mass of the disk is much 
less than that of the planet. We also neglect the effects 
of the PR drag for simplicity; hence our simulations rep¬ 
resent gas-poor planetcsimal or debris disks composed of 
solid bodies 1 m to 100 km across and which are rel¬ 
atively free of dust. The presence of strongly reflect¬ 
ing or emitting dust can markedly affect the appearance 
of a disk, particularly if collision among the dust parti- 
cles is considered . For instance, a study conducted by 
IStark k Kuchnerj (]2009f) shows that the ring structures 
created by the trapping of dust in resonance with a planet 


are smeared out by collisiona l interactions of the dust 
particles. On the other hand, Wyatti (|2QQ5f) has argued 
that the collisional lifetime of dust in debris disk candi¬ 
dates is short enough that dust does not drift very far 
before its destruction. Here we model dynamically cold 
planetcsimal disks which are dust-poor. However, we 
note that the observational characteristics of such disks 
will be dominated by much-smaller dust (which has dif¬ 
ferent dynamics) in dust-rich systems. The mechanisms 
of dust production and removal in these disks are com¬ 
plex and their modeling is outside the scope of this paper. 
Nonetheless, it is reasonable to expect that at least some 
dust-poor systems exist, and it is to those that we turn 
our attention here. 

The simulations are performed at scales appropriate 
to our solar system (i.e., the perturbing planet is placed 
at Jupiter’s semimajor axis, ~ 5.204 AU from the star). 
The physics involved scale with distance, however, and so 
our results are applicable to disks and planets in general, 
even if located at different distances from their parent 
stars. For clarity then, our figures are scaled so that 
the planet is at a unit distance. The exception to the 
scalability of our results is only the timescales for opening 
up the gaps, which are expected to be longer for larger 
systems or ones with less massive central stars. 

3.2. Simulated Debris Disks 

In order to investigate structures in debris disks that 
are caused by MMRs with a planet, we perform simu¬ 
lations of test particles in a flat disk (particle inclina¬ 
tions i' = 0.0°), containing 10,000 particles per 1 AU 
of the disk radial thickness. Running on a single CPU, 
the simulations take 12 hrs to complete for a disk with 
20,000 particles and could last up to 4 days for three 
times the number of particles. The total simulation time 
also depends on the planet’s mass since more particles 
are ejected at the beginning of the simulations with more 
massive planets. The initial particle eccentricities are 
those of the forced eccentricity induced by the planet at 
each particle’s semimajor axis with their apses aligned 
with the planet’s. The forced eccentricity represents the 
eccentricity that the particle orbits are subject to due to 
the simple presence of the perturbing planet. By setting 
the initial conditions of the disk to this value, we create 
a disk which will be minimally perturbed by the planet. 
This choice represents a scenario where the planet and 
the disk have been in the same relative geometry for a sig¬ 
nificant fraction of the age of the system. Other choices 
result in a more heavily perturbed disk, with results that 
are highly dependent on the choice of the initial condi¬ 
tions. Such scenarios might be appropriate to cases of 
recent planetary migration but are not considered here. 
We note that for the case of a planet on a perfectly cir¬ 
cular orbit, the forced eccentricity is also zero and the 
disk particles are started on circular orbits. 

For simplicity, we assume that the planetary system 
contains one planet only and we choose the planet ’s semi¬ 
major axis to be that of Jupiter (a = aj = 5.204 AU), 
though as noted earlier the choice of this scale is arbitrary 
and the results apply equally to disk-planet systems of 
all sizes. Furthermore, in all of our simulations, we as¬ 
sume that the planet orbits in the same plane as the disk 
(i.e. i = 0.0°). However, we also tried some simula¬ 
tions with a small planet inclination (i = ij = 1.304°, 
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Figure 2. Maximum libration widths for 2:1, 3:2, and 3:1 main belt resonances with Jupiter. The vertical dotted lines are the locations 
of nominal resonances (calculated using Equation f2]i and the horizontal dotted line shows Jupiter’s eccentricity. 


Jupiter’s inclination). The results show no noticeable 
difference between the two cases. Therefore, even if the 
planet does not orbit exactly in the disk plane, we still 
expect the same features and to obtain similar results for 
small orbital inclinations. 

The simulations are performed for the case of both ul¬ 
terior and exterior resonances; and in both cases, the disk 
is placed 1 AU away from the planet (2.204 < anisk < 

4.204 for the interior resonance and 6.204 < a Disk < 

12.204 for the exterior case). The values for the disks’ 
inner and outer edges are chosen such that the three reso¬ 
nances being considered, the 2:1, 3:2, and 3:1, fall within 
the debris disks. A range of planet masses is used, going 
from as small as 1.0 M® to 9.0 Mj , where A/® and Mj 
are the mass of the Earth and Jupiter respectively. Fur¬ 
thermore, to study the effect of the planet’s eccentricity 
on MMR gaps, two different planet eccentricity values 
are considered: e = 0.0 and e = 0.0489, the eccentricity 
of Jupiter (ej), though we leave further investigations of 
eccentricity effects to a follow-up paper. 

4. RESULTS 

As noted before, we choose our total simulation time 
to be a million years, which we find to be sufficient for 
the disks to achieve a quasi-steady state. Yet we observe 
gaps forming on much shorter timescales for more mas¬ 
sive planets (e.g., only 100,000 years, for a 5 Jupiter-mass 
planet exterior to the disk). More massive planets open 
gaps more quickly as would be expected, but we have 
not investigated this trend in the present work. At the 
end of the simulations, the disks are examined for struc¬ 
tures, particularly those that would be observationally 
discernible in a telescopic image. Moreover, in order to 
be able to compare the widths of the gaps that are pro¬ 
duced in the simulations with the analytical calculations 
of the maximum libration widths through Equation^ we 


make histogram plots of the number of particles per he¬ 
liocentric distance. The disk is divided azimuthally into 
four equal segments since the MMR gaps turn out to be 
azimuthally asymmetric, particularly in the case of the 
perturbing body having zero eccentricity. 

It must be noted that although we calculate what the 
resonance widths should be analytically for the 3:2 and 
the 3:1 resonances as well as the 2:1, due to the nar¬ 
rower gaps these resonances produce we only measure 
the widths of the 2:1 gaps i n our simulations. However, 
as we shall discuss in Section 15.hi we do observe a narrow 
feature at the 3:1 resonance. 


4.1. Interior Resonance 

Interior resonance refers to the case in which the planet 
is exterior to the disk and hence the resonances occur in¬ 
terior to the planet’s orbit. In this case, the 2:1 MMR 
opens gaps in two regions: one at inferior conjunction 
with the planet and one at opposition. This can be seen 
in Figure |3(a)| which shows the case of a 4.0 M j planet 
on a circular orbit gravitationally interacting with the 
disk. The locations of the three prominent resonances 
are marked; the gap occurs at the 2:1 MMR with the 
planet. If we were to follow the double-arced gap over 
time, we would see it move around the star at the same 
rate as the planet. The appearance of such a gap in an 
image of a disk would both indicate the presence of a 
planet as well as its location (that is, along a fine drawn 
through the deepest parts of each arc), though which side 
of the disk the external perturber is on cannot be unam¬ 
biguously determined. This double-arc shape is a result 
of the removal of particles from the 2:1 interior reso nan ce 
(see for example Fig. 8.4(a) of [Murray &; Dermottl11999ft . 

In order to measure the width of the gap seen in the 
simulation, we make a histogram for each of the four 
colored segments (separated by intervals of 7r/2, start- 
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ing with -7r/4 < 6 < 7r/4) shown in Figure |3(a) for the 
number of particles as a function of thei r dist ance from 
the central star. This is shown in Figure |3(bJ where the 
colors of the first fou r histo grams correspond to the same 
colors used in Figure 3(a) while the last histogram shows 
the distribution of all the particles from the four seg¬ 
ments put together. For simplicity, we will concentrate 
on the gap in the r egion closest to the planet (mar ked i n 
red in Figurc|3(a)|and the top histogram of Figure|3(b)[). 

The dotted lines in Figure 3(b) are the theoretical lo- 
cations of the three resonances that are being considered, 
calculated using Equation 0 while the dashed lines show 
the theoretical width of each resonant gap, obtained from 
Equation [4] In addition to the gaps that appear at the 
2:1 MMR, the histograms also show a slight decrease in 
the number of asteroids in the region corresponding to 
the 3:1 MMR with the planet, although it is considerably 
narrower compared to the 2:1 resonance. The decrease 
at the 3:1 resonance in the sector nearest the planet is a 
difference of about 50% which is 5 times the Poisson er¬ 
ror for the bin and hence is statistically significant at this 
resolution. On the other hand, the 3:2 resonance which 
is closer to the planet cannot be seen since large numbers 
of particles in the outer edge of the disk have been re¬ 
moved due to strong gravitational interactions with the 
planet. To obtain a measure of the width of the 2:1 gap, 
we make a Gaussian fit to the histogram where the gap 
appears for comparison to the analytic results which we 
discuss in Section HT31 

Figure [TTJ belongs to the case in which the planet per¬ 
turbing the disk has zero eccentricity. When the planet’s 
eccentricity is increased to 0.0489, features are seen in 
the disk that are similar to those of the zero eccentricity 
case. Additionally, the feature at the 3:1 MMR broadens. 
This feature, unlike the 2:1 gap, does not orbit the cen¬ 
tral star at the same rate as the planet. This is shown in 
Figure [4]for a planet with e = ej = 0.0489. We attribute 
this to the generation of tightly wound spiral waves at 
a Lindblad resonance in the disk, a phenomenon we ex¬ 
plore in more detail in Section [5.61 

4.2. Exterior Resonance 

Resonance structures are different depending on 
whether the disk is internal or external to the orbit of 
the planet. For the case of exterior resonance (i.e., the 
planet being interior to the disk), we find the shape of 
the gap to be different. Instead of forming two arc-like 
sectors, the gap formed due to resonant interactions with 
a planet on a circular orbit is a single arc whos e cen ter 
aligns with the planet. This is shown in Figure|5(a)|for 
a planet with M = 4.0 Mj. The arc again co-rotates 
with the planet and provides an indicator of the location 
of the perturbing body even if it were not visible in an 
image of the disk. We note again that the most promi¬ 
nent gap appears where the 2:1 MMR with the planet 
is (shown by the filled square). The different shape of 
the gap can also be explained in the same manner as 
for the i nterior resonance a n d is illustrated in Figure 
8.4(c) of IMurrav h Dermottj (|1999 1. Thus the differing 
gap shapes for interior versus exterior resonances allow 
a great deal of information about a perturbing planet to 
be gleaned from images of a disk that displays gaps, even 
if the planet itself remains unseen. 

Following the same analysis that is done for the interior 


resonance, we make a Gaussian fit to the gap that can be 
seen in the histogram for the r egion closest to the planet 
(see top panel of Figure |5(b)| to obtain the width and 
the mean location of the 2:1 MMR gap. We d isc uss t he 
results of these measurements in Sections 15.31 and 15.41 

Figure [5] is obtained when the perturbing body (the 
planet) has no eccentricity. When the planet’s orbital 
eccentricity is increased to that of Jupiter, the arc in 
the 2:1 gap remains easily visible but extends further, 
becoming more annular in shape. This is shown by Fig¬ 
ure 03 In addition, the 3:1 feature observed in the case 
of the interior resonance (Figure 0 also becomes more 
prominent, and is explained in Section [5.61 

5. DISCUSSION 

The main purpose of this study is to provide the means 
of determining the properties of a perturbing planet in 
the case where an MMR gap is observed in a debris disk. 
This includes both the characterization of detected plan¬ 
ets, as well as providing information on the mass and 
location of planets which remain as-yet unseen. We con¬ 
sider now how the results of our simulations provide in¬ 
sight into this issue. 

5.1. Asymmetries in MMR Gaps 

Our simulations suggest that the shape of an MMR gap 
is different for interior and exterior resonances. There¬ 
fore, by looking at the shape of a gap in a planetesimal 
disk we are not only able to immediately determine which 
side of the disk the perturbing planet is, the azimuthal 
asymmetry of the MMR gaps also allows us to easily 
distinguish between gaps formed due to MMRs with a 
planet that lies outside the gaps and those formed by 
planets in the gaps for which the gaps arc azimuthally 
symmetric (such as the ones seen in the HL Tau disk). 
The only exception, according to our simulations, is the 
case of exterior resonance with a planet having non-zero 
orbital eccentricity for which we obtain nearly annular 
gaps at the 2:1 resonances. We note here that it is possi¬ 
ble that increasing the eccentricity of the planet would re¬ 
sult in the gap to eventually become completely annular, 
although we do not see that in a single sample simulation 
with e = 0.1. However, even in this case, we still observe 
an asymmetric feature at the 3:1 MMR with the planet. 
The presence of this additional gap at the 3:1 resonance, 
which occurs for both interior and exterior resonances, 
is indicative of the planet having non-zero eccentricity. 
We leave examining the effect of planet eccentricity on 
debris disk struct me to a future paper. 

5.2. Minimum and Maximum Detectable Planet Masses 

We find that planets as small as M = 1.0 M® can pro¬ 
duce MMR gaps (see Figure0, but the practical lower 
limit on the mass of the planet that can open a detectable 
gap depends on the resolving capabilities of the observa¬ 
tional facility taking the images. We are not aware of any 
disk with features meeting the above descriptions (i.e., 
disks with observed azimuthally asymmetric MMR gaps) 
that have yet been reported, but given the resolutions ob¬ 
tained by current facilities, we expect such features to be 
discovered in the near-future. For ill ustration, th e disk 
around HL Tau has a radius of 80 AU (|Kwon et al.||2011|) 
and the recent ALMA image of the disk shows gap fea¬ 
tures as small as 5 AU across (|Tamayo et alj|2015() . For 
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(a) Interior resonances result in formation of two (b) Distribution of particles in the disk and MMR structures due to a planet exteri or to 
arcs of gaps in the disk. The dashed line shows the disk. The colors in the first four panels correspond to the same colors in Figure[3(aj| 
the periastron of the planet's orbit and the sym- for the different segments (from top to bottom: regions separated by intervals of tt/ 2, 
bols represent the three resonances considered starting with —rr/4 < 0 < rr/4). The last panel shows the overall nmnber distribution 
while the filled circle is the planet. Different col- of the particles in the disk. The dotted lines are the theoretical locations of the 2:1, 
ors represent the different segments for which we 3:2, and 3:1 resonances with the dashed lines defining the width of each gap calculated 
make histograms of th e dist ributions of particles analytically. A Gaussian fit is made to the top histogram where the gap is to obtain a 
in the disk (see Figure [3(b) | |. measure of the gap width from the simulations. The bin size is 0.006 a, in the unit of the 

planet’s orbital radius. 


Figure 3. Simulation showing interior resonance structures formed by interactions of planetesimals in the disk with a 4.0 Mj planet on 
a circular orbit exterior to the disk. 


comparison, a 3 Jupiter-mass planet on a circular orbit 
just outside the edge of a planetesimal disk of similar 
size as the HL Tau disk would create a 2:1 resonance gap 

that is W = (0.013 x 3.0 + 0.020) x 80.0 ~ 4.7 AU wide, 

comparable to the observed gaps that have been associ¬ 
ated with planet formation, and easily distinguishable. 
It must be noted that ALMA probes particles that are 
~ 1 mm and smaller in size, but here we model much 
larger bodies. Nonetheless, millimeter-size particles are 
nearly unaffected by the PR. drag and thus we obtain 
similar results when we consider disks that are entirely 
composed of 1 mm dust and when colhsional dust pro¬ 
duction is ignored. Yet future work is needed to study 
MMR gaps when collision among millimeter and sub- 
millimeter particles and radiative forces are taken into 
consideration. On the other hand, if the planet becomes 
very massive, the edge of the disk thins under its pertur¬ 
bations as larger mass planets excite the particles and 
move them to higher eccentricity orbits. In our simula¬ 
tions, the disk becomes too heavily eroded to distinguish 
the 2:1 MMR at 6 - 7 Mj , making this the practical up¬ 
per limit for the use of this technique for characterizing 
extrasolar disk-planet systems. 


5.3. Libration Width vs. Planet Mass 

More massive planets are expected to open wider MMR 
gaps in disks. In Figure[S]we illustrate this effect by plot¬ 
ting libration width versus planet mass for interior (left) 
and exterior (right) resonances with planet’s eccentric¬ 
ity e = 0.0, shown by the top two panels, while the 
bottom two panels show our results when the planet’s 
eccentricity is increased to e = 0.0489. The different col¬ 
ors and symbols used and their least-square linear fits in 
each panel represent results obtained analytically using 
Equation [jl] and from the simulations. It is clear that 
particularly in the case of interior resonances (left two 
panels), there is nearly perfect agreement between the 
2:1 MMR gap widths obtained from the simulations and 
the ones calculated analytically while there is an offset 
between the measured and analytic results in the case 
of exterior resonances (right panels). Furthermore, as 
expected, increasing the planet’s mass also increases the 
width of the MMR gap in all cases. Thus the mass of 
the perturbing planet, seen or unseen, can be estimated 
from the width of the resonant gap it produces. 

Equations [5] through 0 are obtained from least square 
fits to the values we get from our simulations (sec 
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(a) Same as Figure [3(a) | except that the planet’s (b) Histogram plots for the disk shown in Figure |4(a)| also indicate the presence of an 

orbital eccentricity is increased to ~ 0.05. The extra gap at the 3:1 resonance with the planet. 

extra feature at the 3:1 MMR is likely launched 

by Lindblad resonances and appears whenever 

the planet is given non-zero eccentricity. 


Figure 4. Same as Figure[51 but with planet’s eccentricity increased to e = ej = 0.0489. 


Figure [8j. These equations allow one to easily calculate 
the planet’s mass if the 2:1 MMR gap width can be 
measured observationally. To simplify the calculation, 
we normalize our results again, this time making the 
observed mid-location of the gap equal to one unit of 
distance. We do this so that the equations that are 
provided from this point forward can also be applied to 
cases where the planet is unseen but an MMR gap is 
detected in the debris disk. 

(A) Interior Resonance with e = 0.0: 

M= 0l5 (W ''"°- 012)l (6) 

(B) Interior Resonance with e = 0.0489: 

M= om {w °~ Q - 025) ’ (7> 

(C) Exterior Resonance with e = 0.0: 

M = o^9<W- - 0.009), (8) 

(D) Exterior Resonance with e = 0.0489: 

e 08 ( ^-°'° 20) 


where M is the planet mass (in Jupiter masses, Mj) 
and W 0 is the observed width of the gap (in units of the 
distance between the star and the observed gap). 

The error bars are calculated by taking three main 
sources of uncertainty into account that are added in 
quadrature: 

1) The Gaussian fit to the histogram is made by least- 
square fitting using Interactive Data Language (IDL) and 
the goodness of fit is recorded as one source of uncer¬ 
tainty. This value is generally small in the examples we 
tried. 

2) Since the Gaussian fit is made to points that mark 
half the bin size, the measurements have uncertainties 
that are affected by the choice of the bin size (taken 
to be 0.006 times the planet’s distance from the star in 
this work). However, the uncertainty in each bin position 
goes as \fn , where n is the number of particles in the bin; 
and since we chose the bin size such that on average each 
bin contains not more than 0.25% of the total number of 
particles in the disk, the uncertainty due to the finite size 
of the bin is usually not significant. 

3) Examining the shapes of the MMR gaps in our simu¬ 
lations revealed that gaps are not always perfectly Gaus¬ 
sian in shape. Thi s is i llustrated, for instance, in the 
top panel of Figure |3(b)| where the gap is higher at one 
end. Therefore, in order to find the width of the gap, 
the Gaussian fit is made three times by normalizing to 


M = 


( 9 ) 
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3:2 2:1 3:1 



* <•) Heliocentric Distance (a) 


(a) Exterior resonance with a planet on a circu- (b) Distribution of particles in the disk for each segment marked in Figure [5(a)| 
lar orbit results in the formation of a gap that 
appears as a single arc in the region of the disk 
closest to the planet. 


Figure 5. Structures formed by resonance interactions of planetesimals with a 4.0 Mj planet on a circular orbit interior to the disk. 


either side and also without normalization. The stan¬ 
dard deviation between the three values obtained is then 
taken as the uncertainty in the gap width and is the dom¬ 
inant source of uncertainty in our calculations. It must 
be noted that our results for the gap widths are shown in 
terms of each gap’s standard deviation while an observer 
measuring the width of the gap might define the edges by 
measuring the peak brightness in the disk near the edge, 
then locating the radius in the gap at which the bright¬ 
ness of the disk is half that value. This met hod of us- 

S ' ’ 2 “half-maximum radius” is used by iChiang et al.l 
for simulations of the Fomalhaut disk. Neverthe¬ 
less, the two quantities are related by a simple formula 
showm by Equation 1101 

FWHM = 2V2 x ln2 a, (10) 

where FWHM is the full width at half maximum and o 
is the standard deviation ( W Q in our equations). 

5.4. Resonant Location Versus Planet Mass 

Examining the 2:1 MMR gaps in our histograms sug¬ 
gest that there is a shift in the position of the gap (the 
peak of the Gaussian fit) from the nominal resonance 
location (the dotted fines on the histograms). In fact, 
our simulations show r that the offset from the gap’s theo¬ 
retical position is proportional to the mass of the planet. 
This is showm in Figure®where the top and bottom pan¬ 
els are for interior and exterior resonances, respectively. 


The red squares show the location of the 2:1 resonance 
with a zero-eccentricity planet while the diamond sym¬ 
bols in green are for the case of higher eccentricity planet 
with e = 0.0489 and the dotted fines denote the location 
of the nominal resonance calculated using Equation [2] for 
interior and exterior resonances. The uncertainties in the 
observed gap locations are calculated in the same manner 
as those for the gap widths, explained in Section 15. ill 
In addition to a shift in the resonance location, we note 
from Figure® that the MMR gaps always tend to shift to¬ 
ward the planet as can be seen from the negative slope in 
the top panel and the positive slope in the bottom panel. 
We find that there is small contribution from the planet’s 
eccentricity in shifting the locations of the MMR gaps. 
In addition, our results indicate that the location of the 
2:1 gap agrees better with the theoretical prediction for 
the higher eccentricity planet in the exterior resonance 
case. 

5.5. Finding the Planet's Semimajor Axis based on the 
Observed Gap Width and Location 
If the mid-location and the width of an MMR gap can 
be obtained through observations, this information can 
be used to calculate the planet’s semimajor axis. This 
could be useful for calculating the orbit of a knowm planet 
or determining the location of an unseen one. This cal¬ 
culation is done by first measuring the distance from the 
star to the center of the gap, a! 0 and then finding the gap 
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(a) Same as Figure |5(a)| except that the planet is 
given some eccentricity (e ~ 0.05). Similar to the 
case of the interior resonance with the pertur bing 
planet on a non-circular orbit (Figure [4(a)) , an 
extra gap appears in the disk with location corre¬ 
sponding to the 3:1 MMR with a planet interior 
to the disk. 



Heliocentric Distance (a) 

(b) Histograms corresponding to the different segments of the disk shown in Figure [6(a)| 


Figure 6. Same as Figure[51 but with planet’s eccentricity increased to e = ej = 0.0489. 


width in units of the observed star-gap separation. Then 
Equations [g] through 02 can be used to obtain the planet ’s 
mass. The theoretical location of the gap, a', (i.e. the 
nominal resonant location) can then be found through 
the following equations: 

(A) Interior Resonance with e = 0.0: 


a' = -0.003 M + 0.002 + a' Q , 

(ii) 

(B) Interior Resonance with e = 0.0489: 


a! = -0.002 M - 0.015 + a! 0 , 

(12) 

(C) Exterior Resonance with e = 0.0: 


a' = 0.004 M - 0.001 + a! 0 , 

(13) 

(D) Exterior Resonance with e = 0.0489: 


a = 0.001 M - 0.002 + a oy 

(14) 

where M and a 1 have units of Jupiter mass ( Mj ) and 
a' Q , respectively. Finally, the semimajor axis of the 
planet can be calculated using Equation [2J Equations 11 11 
through [T4] were obtained by subtracting the observed 
gap location from its theoretical location in our simula¬ 
tions. This is shown by Figure [TT71 which can be used to 


obtain the theoretical MMR gap location from which the 
planet’s semimajor axis is found. 

5.6. The Effect of Planet's Orbital Eccentricity and 
Lindblad Resonances 

When the planet is given a small but non-zero ec¬ 
centricity, the 2:1 MMR properties remain largely unaf¬ 
fected. However, tightly wound spiral waves originating 
from the 3:1 MMR appear in some cases. These features 
are potentially valuable sources of information about the 
disk’s properties, but are far more challenging to detect 
in real telescopic images than the gaps associated with 
the MMRs that we have been discussing so far. In fact, 
they often appear at scales below those shown in the pre¬ 
vious figures, their primary visible tracer being a narrow 
arc-like gap at the 3:1 interior or exterior resonance. 

We interpret these features as force d eccentric ity waves 
originating at Lindblad resonances (iShul 1 19841) . These 
waves are similar to those seen in Saturn’s rings (e.g. 
[Holberg et al.|ll982j; ILane et apl982) and have the char¬ 
acteristic decre ase in wavelength as one m oves away from 
the resonance ([Murray & DcrmottJ 119990 . This is illus¬ 
trated by Figure [TTI which shows a series of density waves 
originating from the 3:1 resonance. 

These waves are seen to be produced primarily at the 
3:1 MMR for both internal and external planetary per- 
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(a) No MMR gap is obvious in this simulated 
disk due to the small mass of the perturbing 
planet (M = 1.0 A/®) 


3:1 2:1 3:2 



Heliocentric Distance (a) 

(b) Histograms for the four segments are noisy, but putting them together averages the 
noise out and reveals two narrow gaps at the 2:1 and 3:2 MMR with the planet. 


Figure 7. Though MMR gaps are not obvious in this example with a 1.0 A/® perturbing planet, examining the last histogram reveals 
two dips at the 2:1 and 3:2 MMR with the planet. Such small gaps might not be detectable in a telescopic image of the disk; nevertheless 
our simulations show that even an Earth-mass planet is able to open gaps in disks through MMRs. 


Interior Resonance: Libration Width Vs. PUnci Mass *1101 » 0.0 




Exterior Resonance: Libration Width Vs. Planet Mass wHi e^ » 0.0 



Esteriar Resonance: Libration Width Vs. PUnel Mass with 0.0489 



Figure 8. Libration width versus planet mass. Left panels: interior resonance with e = 0.0 (top) and e = 0.0489 (bottom). Right 
panels: exterior resonance with e = 0.0 (top) and e = 0.0489 (bottom). W and W 0 are the theoretical and the simulated width of the gap, 
respectively, and are given in units of the observed mid-location of the gap (aj,). The vertical axis is the full width of the gap, 2cr of the 
fitted Gaussian for the simulations. 
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Interior Resonance: Resonance Location Vs. Planet Mass 



Planet Mass.M (My) 


Exterior Resonance: Resonance Location Vs. Planet Mass 



Planet Mass.M (M ; ) 


Figure 9. The effect of planet's mass and eccentricity on the resonant location for interior (top) and exterior (bottom) resonances. The 
red squares and the green diamonds show the theoretical gap locations for cases with e = 0.0 and e = 0.0489, respectively while the dotted 
lines show the locations of the observed gaps. 


Intorior Resonance: Gap Shift Vs. Planet Mass 
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Figure 10. The shift in the gap location versus the planet's mass for interior (top) and exterior (bottom) resonances for two different 
planet eccentricities. The symbols used are the same as those in Figure0 a' and a' Q are the theoretical and observed location of the gap, 
respectively. 
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turbers and are likely the reason we see an extra fea- 
ture at the 3:1 resonance with the planet whenever the 
planet is given non-zero eccentricity in our simulations 
(see Figures [I] and [t>]). These coincide with the location 
of the m = 2, k = ±1, p = 0 inner/outcr Lindblad 
resonances, which arc associated with similar waves in 
Saturn’s rings. Though such structures are a rich source 
of information about the planet and the disk itself, the 
propagation of waves in real syste ms depends on effects 
such as self-gravity and collisions flFridmanl 11999f) which 
are not modeled here, and we will leave their examination 
for future work. 

Although the 3:1 MMR gaps are narrower and more 
difficult to measure than the 2:1 gaps, if they can be de¬ 
tected, their relative location with respect to that of the 
2:1 resonance can be used to also distinguish interior from 
exterior resonances as the 3:1 resonance gap is formed 
farther from the planet than the 2:1. We will leave de¬ 
tailed examination of resonant interaction between disk 
particles and higher eccentricity planets to a follow-up 
paper. 


5.7. Disk Optical Depth and Gap Contrast 

Though the gap widths can easily be measured in sim¬ 
ulations, this may prove more difficult observationally. 
In particular, the optical depth in the gap versus the 
disk as a whole will determine the amount of contrast 
in the image. Figure \TZ\ shows the ratio of the average 
disk surface density to that at the deepest part of the 2:1 
MMR gap (the “contrast”), as a function of planet mass. 
Note that the contrast we refer to here is the contrast in 
planetcsimal surface density, not dust surface density, be¬ 
cause we are considering dust-poor systems. Because of 
the difficulty of defining the edges of the gaps, we do not 
try to define an edge-to-center contrast. Though higher 
mass planets form larger gaps, the contrast is somewhat 
diminished by particles “bleeding” in from the edges of 
the gaps. There is a trend for the contrast to increase 
with mass for a planet that is interior to the disk (Figure 
1121 bottom panel) but there is little effect for an external 
perturber (Figure[12] top panel). Thus the gap contrast 
itself can provide a measure of the planet mass in some 
cases, though this approach is observationally more dif¬ 
ficult. Nonetheless, the contrast remains at large values 
through the range of planet masses considered here, indi¬ 
cating that the MMR features discussed here should be 
detectable if the disks themselves are. 

We note from our results that the contrast grows al¬ 
most linearly with planet mass, except for exterior res¬ 
onance with a planet on a circular orbit in which case 
there seems to be an exponential trend. However, we 
feel it is unwarranted to fit an exponential to this case 
for two reasons: First, the appearance of exponential 
growth is only present for one of the four cases (exterior 
resonance with a planet having e = 0) and rests on only 
the two rightmost points. Moreover, those two points 
are near the largest masses beyond which the gaps dis¬ 
appear due to heavy erosion of the disk edge; and so the 
graphs should not be extrapolated beyond the presented 
maximum value. 

Simu lated ob servat ions of the disks shown ear lier in 
Figures|3(a)|and|5(a)[ are illustrated in Figures [T31and[Til 
assuming an inverse-square dependence of particle emis¬ 
sion on distance from the central star. It must be noted 


again that in this study, we assume that disks are opti¬ 
cally thin and are largely free of collisionally produced 
or other sources of dust. These simple figures are free of 
additional noise that would certainly be present in real 
observations. The pixel size is chosen to match ALMA’s 
highest resolution at 1 mm wavelengths (15 mas for the 
16 km configuration) at the distance of HL Tau (140 pc) 
if the disk’s outer radius is 100 AU. Since the MMR gaps 
can still be seen in these images, we conclude that the 
MMR features described here are nominally within the 
reach of current observational facilities. 

6. SUMMARY AND CONCLUSIONS 

We investigated the dynamical effects of a planet on a 
planetcsimal disk through MMRs for both interior and 
exterior resonances. Our purpose is to use the observed 
properties of MMR structures to characterize the planet 
producing them, even if that planet remains as-yet un¬ 
detected. 

Structures arising from MMRs can be highly diagnos¬ 
tic of the properties of the planet disturbing the disk 
particles. MMR gaps become wider as the planet’s mass 
increases; there is a linear relationship between the gap 
width and the perturbing planet’s mass that agreed well 
with theoretical calculations. Therefore, measurement of 
the width of a 2:1 MMR gap would help determine the 
perturbing planet’s mass, even if it remains unseen. We 
find gaps at the 2:1 and 3:2 MMRs even for a planet as 
small as 1.0 A/®, although their small widths make them 
observationally more challenging to detect than those at 
Jupiter masses. On the other hand, at planet masses be¬ 
yond 6 — 7 Mj, the resonance structures are destroyed as 
the disk is eroded by the planet ’s growing Hill Sphere. 

We found an offset in the gap’s position in the disk 
from the nominal resonant location with more massive 
planets causing a larger shift in the observed location of 
the gap. Thus if the planet’s location is already known, 
the shift from the theoretical location of the 2:1 MMR 
gap can be used to confirm the planet’s mass, which can 
alternatively be calculated using the gap width. On the 
other hand, for cases in which the planet remains un¬ 
detected, we proposed a set of equations that take the 
planet’s mass, calculated using the gap width, to deter¬ 
mine the planet’s location based on that of the observed 
2:1 gap. 

We further extended our studies to simulations of sys¬ 
tems in which the perturbing planet has non-zero orbital 
eccentricity. In this case, disk structures due to MMRs 
become more complicated and the 2:1 MMR gaps formed 
by an internal perturber become more annular in shape. 
Moreover, an extra arc-like feature was seen at the 3:1 
MMR with the planet when the planet’s eccentricity was 
increased and is associated with spiral waves generated 
at Lindblad resonances. Thus the appearance of an arc¬ 
shaped gap at the 3:1 resonance with the planet is in¬ 
dicative of the planet having non-zero eccentricity. If 
detected, it could also be used to distinguish interior 
from exterior resonance. This would be helpful, partic¬ 
ularly for the case of exterior resonance with a non-zero 
eccentricity planet for which we found the 2:1 gaps to 
be azimuthally symmetric. We obtained nearly perfect 
agreement between the 2:1 gap widths measured analyt¬ 
ically and through our simulations for a planet exterior 
to the disk with 0 or 0.0489 orbital eccentricity. On the 
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Figure 11. Spiral patterns appear commonly in our simulations when the perturbing planet has non-zero orbital eccentricity. These are 
likely forced eccentricity waves originating at Lindblad resonances. Our simulations indicate that these waves are generated at a location 
corresponding to the the 3:1 MMR with the planet, marked by a diamond on this figure. 
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Figure 12. MMR gap contrast vs. planet’s mass (in Jupiter mass) for interior (top) and exterior (bottom) resonances. 
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Figure 13. Same as Figure |3(a) | for interior resonance except that particles in x and y are binned and assigned a color based on the total 
emission from each bin, normalized to the peak intensity. Darker colors correspond to less emission. 
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Figure 14. Same as Figure IT51 but for the case of exterior resonance. The same disk is shown in Figure [5(a)| 
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other hand, the difference between the two measurements 
seems to grow with mass for a planet on a circular orbit 
interior to the disk while a systematic shift is seen when 
the planet’s eccentricity is increased to that of Jupiter’s 
in this case. All the simulations we report on here were 
performed for flat disks with the perturbing planet in the 
same orbital plane, but we found no significant difference 
in our results when the planet, and thus the particles, 
were given a small orbital inclination. 

The results of our simulations indicate that the shapes 
of the gaps opened by the 2:1 MMRs are different for in¬ 
terior versus exterior resonances, with the former making 
two sectors at the planet’s (inferior) conjunction and op¬ 
position while the latter forms a single arc at (superior) 
conjunction. Since direct detection of extrasolar plan¬ 
ets still remains observationally challenging, detection of 
such structures in a planetesimal disk allows one to not 
only infer the presence of an unseen planet, the two dis¬ 
tinct gap shapes would also make it possible to easily de¬ 
termine the relative location of the planet with respect to 
the disk and to distinguish MMR gaps from azimuthally 
symmetric gaps formed by a planet that is embedded in 
the disk. As ALMA and other facilities continue to ad¬ 
vance the frontiers of extrasolar planetary science, the 
ability to detect and characterize unseen planets based 
on their effects on a more-easily observable disk will be¬ 
come an increasingly powerful tool. 
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